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Abstract: In this work we present a minimal parametrization of the light-cone distribu¬ 
tion amplitudes of the baryon octet including higher twist contributions. Simultaneously 
we obtain the quark mass dependence of the amplitudes at leading one-loop accuracy by the 
use of three-flavor baryon chiral perturbation theory (BChPT), which automatically yields 
model-independent results for the leading SU(3) flavor breaking effects. For that purpose 
we have constructed the nonlocal light-cone three-quark operators in terms of baryon octet 
and meson fields and have carried out a next-to-leading order BChPT calculation. We were 
able to find a minimal set of distribution amplitudes (DAs) that do not mix under chiral 
extrapolation towards the physical point and naturally embed the A baryon. Additionally 
they are chosen in such a way that all DAs of a certain symmetry class have a similar quark 
mass dependence (independent of the twist of the corresponding amplitude), which allows 
for a compact presentation. The results are well-suited for the extrapolation of lattice data 
and for model building. 
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1 Introduction 

Due to the unstable nature of the weakly decaying hyperons there are no scattering experi¬ 
ments with hyperons in the initial state. However, they naturally occur in the final state, for 
instance in baryon-antibaryon pair production via electron-positron annihilation e + e^ —> 
BB, in deeply virtual exclusive meson electroproduction 7 *p —> K + A, /v + E°, A'°E + , 
and in decays of heavy quarkonia to baryon-antibaryon pairs like J/\ I/, T —» BB. The 
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standard way to parametrize the nonperturbative information contained in such exclusive 
processes are (transition) generalized parton distributions or ordinary form factors. At high 
momentum transfer the contributions from Fock states containing more than the minimal 
number of partons are power-suppressed and the process can be approximated by a con¬ 
volution of the involved distribution amplitudes (DAs) with the process-dependent hard 
scattering kernel. The requirement of large momentum transfer, the instability of the final 
state hadrons and the fact that distribution amplitudes only occur in convolutions require 
high luminosity and high granularity detectors to extract information on the hyperon DAs 
from experiment. 

Another type of process where hyperon DAs are involved are the exclusive rare decays 
of 6 -baryons, like E*,, A&, E& and finto octet baryons (plus 7 , l + l~, ...). Due to the large 
mass difference one can hope that higher order Fock states are sufficiently suppressed to 
allow for a description by three-quark DAs. Since the bottom baryons are produced with 
increasing rates at LHC and at B-factories worldwide, we have to expect that ever more 
precise experimental results will be available in future, even for rare decays containing 
flavor-changing neutral currents, which are sensitive to new physics. Notwithstanding the 
fact that 6 -baryons are produced at much lower rates than b- mesons, they are not less 
interesting since they allow for an examination of the helicity structure of the b —> s 
transition and thus complement the measurements in the meson sector [ 1 ], As shown in 
refs. [2, 3] there are possible scenarios where deviations from the standard model are not 
seen in the branching ratio of A& —> A l + l~ but only in the A baryon polarization. It 
is therefore mandatory to establish a theoretical basis for the description of such decays, 
and the knowledge of hyperon DAs is one important ingredient. Even the higher twist 
components can yield relevant contributions [4]. Note that constraining the shape of wave 
functions by calculating the moments of the DAs with lattice QCD plays an even more 
important role for hyperons than for nucleons, since experimental bounds are less strict 
than in the nucleon sector. 

A first parametrization of the leading twist contributions in hyperon wave functions was 
already presented in ref. [5]. A complete parametrization (including all contributions from 
higher twist) of baryon-to-vacuum matrix elements was first performed for the case of the 
nucleon in ref. [ 6 ], where it turned out that higher twist contributions can yield substantial 
effects in the baryon sector, since the corresponding normalization constants Af and A 2 are 
large compared to the leading twist wave function normalization constant f N . The same 
procedure has later on been reused in refs. [7, 8 ] to give similar parametrizations for matrix 
elements of the hyperons in the baryon octet, namely E*, E°, S _ , E° and A. Our work 
unifies these different approaches and we find relations between the distribution amplitudes 
for different baryons even if SU(3)/ symmetry is broken. The obtained relations are exact 
including terms up to first order in the quark masses. In this sense we call our results 
model-independent. However, one should keep in mind that higher order contributions 
which lie beyond the accuracy of our analysis are model-dependent indeed, since they are 
affected by the neglection of higher order terms during operator construction and by the 
choice of the regularization scheme. 

As shown in refs. [7, 8 ] one has to introduce six additional DAs if one extends the 
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formalism from the nucleon doublet to the complete baryon octet. Our results show that 
these additional DAs are determined by the eight independent DAs already known from 
the nucleon sector. I.e., if one knows the eight standard DAs (and their dependence on the 
mass splitting between light and strange quarks) for the A and for at least two types of octet 
baryons with nonzero isospin, one can predict all the rest. Using the parametrization given 
in refs. [9-11], where contributions of Wandzura-Wilczek type [12] are taken into account, 
and applying the approximation advocated in ref. [9], where contributions that can mix 
with four-particle operators are systematically neglected, we need only 43 parameters to 
describe the complete set of baryon octet DAs, including their dependence on the splitting 
between light and strange quark mass. For details see section 5.3. This amounts to a 
significant reduction of parameters compared to an ad hoc linear extrapolation without 
the knowledge of SU(3)/- symmetry breaking, which would require 72 parameters for the 
given setup. Therefore our results are useful for the extrapolation of lattice data. In a first 
step it can be checked whether the nontrivial relations between the different DAs that we 
have obtained are realised in lattice simulations. If this is the case to a satisfactory degree, 
one can perform a simultaneous fit to all DAs, which, owing to the significant reduction 
of parameters mentioned above, has much higher precision. Note that the parameters 
occurring in the approximation described above are determined by the zeroth, first and 
second moments of the leading twist DAs and by the zeroth and first moments of twist 4 
DAs, which are, apart from the first moments of the higher twist amplitudes, within reach 
of state of the art lattice simulations (see ref. [13]). 

Let us note that SU(3)j breaking effects can be of considerable strength: e.g. sum rule 
estimations of the symmetry breaking for leading twist wave function normalization con¬ 
stants range from ~ 10%, see ref. [5], to ~ 50%. The latter is obtained if one takes the values 
for / s and f~ from refs. [7, 8]. In ref. [5] it is stated that the impact on the shape of the 
wave function is even larger (at a scale of 1 GeV). One therefore has to expect substantial 
effects also at intermediate scales which are relevant for phenomenological computations. 
We expect lattice QCD calculations to provide quantitative results for SU(3)/ breaking 
effects at physical mean quark mass in the near future. These will provide the means to 
determine the strength of the symmetry breaking terms in our formalism, allowing us to 
draw more definite conclusions. 

At this point we want to highlight a conclusion that can be drawn from our results, 
which is of conceptual importance and also affects the nucleon sector: we find that the 
nonanalytic chiral behaviour of moments of DAs does not depend on the twist of the 
amplitude. Instead, the leading chiral logarithms in the chiral-odd sector are determined by 
the type of amplitude to which the corresponding moment contributes. The ones occurring 
in <& B t (& B i) amplitudes, which will be defined in eq. (5.7), have the same chiral logarithms 
as f B (Af). The odd moments of the leading twist DA therefore behave like Xf instead of 
f B , which is quite contrary to the intuitional expectation. The shape parameters occurring 
in ref. [9] can all be assigned uniquely to one of the two classes, which means that the 
destinction between moments described above is to some extent already present in currently 
used parametrizations. 

This work is organized as follows: In section 2 we present some fundamental definitions 
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to lay the base for the parametrization of the nonlocal three-quark operators in terms of 
baryon and meson fields, which is performed in section 3. A sketch of the leading one-loop 
baryon chiral perturbation theory (BChPT) calculation is given in section 4, where we 
also explain how we have matched our results to the standard DAs given in ref. [6]. In 
section 5 we present our main results. We provide a definition for DAs that do not mix 
under chiral extrapolation and naturally embed the A baryon. The result section is to the 
most part self-contained such that the reader can skip the details of the derivation at will. 
We summarize in section 6. 


2 Fundamental definitions 


There exist various possible realizations of chiral symmetry, which all lead to equal results. 
In the following we only present the definitions we use in this work. For a detailed treatment 
of the effective field theory framework we refer to [14-19]. The pseudoscalar fields are 
contained in 


u 


= ' /z ' =exp (^b A “ <# '“) =exp (^)' 


( 2 . 1 ) 


where A 1 , ..., A 8 are Gell-Mann matrices and Fq is the pion decay constant in the three- 
flavor chiral limit, which corresponds to the convention where F n = Fq + O(. rri 2 K , rn ^) » 


92 MeV. The matrix (f> can be written in terms of meson fields 

7T + K + 

- 77 ^°+ 77 ^ K° 
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^ V2 n ° + V6 1 
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The 3x3 matrix B contains the baryon octet: 


B = 


-7=S 0 + -^=A 

%/2 x /6 

S' 

V “ 


-4=s° + -4=A 

\/2 \/E 

"0 


p 

n 

_Aa 

Ve / 

+ 


(2.3) 


= K p p + K n Tl + Ks-S + 7C s oS U + Kv+S + + +K=o£ U + KaA, 


where the second line defines the matrices kb■ Let us from now on use X e {L,R} and, 
as a convenient notation, L = R and R = L. Where they are not used as an index, L and 
R are meant to be elements of SU(3 )a/_r- Defining ur - u and ul = u ^ the transformation 
properties of meson and baryon fields under chiral rotations read 


u x Xu x K ] = Ku x X ] , (2.4a) 

B KBK f , (2.4b) 


with the so-called compensator field K, which is a common, nonlinear realization of chiral 
symmetry [17, 20]. The covariant derivative acting on a baryon field is defined as 

D li B = d li B + [T li ,B], (2.5) 
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Vr 

Vr 

Vr 

Vr 

1 

1 

-1 

1 

1 

75 

-1 

-1 

-1 

1 

7/7 

(-I)m 

1 

1 

-1 

7/t75 

-(-IV 

-1 

1 

-1 


(- 1 U- 1 ), 

1 

1 

1 


Table 1 . The constants q j?, , rfe and 77 ^ characterizing the symmetry properties of the elements 

of the Clifford algebra, where (-1) M is 1 for q = 0 and -1 for /x = 1,2,3. 


where T^ is called the chiral connection and is given by 

T/j. = ^{u^d^u + ud^v)) . (2.6) 

The chiral vielbein u tl and the quark mass insertions X ± are defined as 

U/j, = ifu'd^u - ud pvt) , (2.7a) 

X± = ± ux’u , (2.7b) 


where x = includes the quark mass matrix, and transform under chiral rotations as 

follows: 


u,j, Ku^K^ , ( 2 . 8 a) 

X± K X J< ] ■ (2.8b) 


Finally we define for the elements of the Clifford algebra in a unitary representation 


r = Vrlo^lo , 

(2.9a) 

r T = ? 7 p CTC , 

(2.9b) 

r f = ?7r7 0 T7o , 

(2.9c) 

10 

[-H 

Win 

11 

t-H 

(2.9d) 


where C = i'y 2 'y° is the charge conjugation matrix. The different t/s are collected in table 1. 

3 Operator construction 

In this section we will construct the light-cone (n is a lightlike four-vector) three-quark 
operator 

qa( a in)q b p(a 2 n)q*(a 3 n ) , (3.1) 

in terms of baryon octet and meson octet fields. The antisymmetrization in color indices 
(which makes the operator a color singlet) and the Wilson lines connecting the quark fields 
(providing gauge invariance) are not written out explicitly, a, b, c are flavor and a, j3, 7 
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Dirac indices. Note that there are many possible parametrizations owing to the freedom 
of choice one has by neglecting higher order effects. The task is therefore not only to find 
a parametrization, but to find one that is most convenient for the loop calculation to be 
performed and can be easily matched to the standard decomposition given in ref. [6]. For 
the parametrization of the nonlocal operator one needs functions, where the moments of the 
functions play the role of low energy constants (LECs). For the parametrization presented 
below these functions can be easily matched to standard distribution amplitudes. 

3.1 Symmetry properties 

To perform the construction of an operator within the effective theory we have to know its 
symmetry properties. To make use of chiral symmetry it is convenient to split the quark 
fields in left- and right-handed parts 

ql(ain)q b p(a 2 n)q°(a 3 n) = O a ^ olM {a l ,a 2 ,a 3 ) + Of£ a p 1 (a 1 ,a 2 ,a 3 ) 

+ ORL,aP 7 ( a l> 02,03) + ®Lli,aP'y(ai,a 2 ,a 3 ) 

+ ^RL,'yap ( a 3^ a l ^ 2 ) + C , L C R,7a/3( a 3i a l)°2) 

+ £ > RL,0'ya( a 2> a 3,Ol) + ®*LR,/3ja(. a 2, a 3, a i) , 

where the operators Oxy for 1,7 e {L,R} are given by 

®XY,aP'y( a l> a 2> a 3) = q > X,a( a 3 n ) ( lX,p( a 2 n ) ( lY,'y( a 3 n ) ; 

where the left-/right-handed quark fields are defined as qL/R = 7 l/rQ with the projection 
matrices 7 l/_r = (1 t 7s)/2. These operators can be characterized by their transformation 
properties under parity transformation (p), charge (c) and hermitian (f) conjugation and 
chiral rotations (y): 

U P b 

®XY,a/3"/( a l’ a 2 i a 3) * (7o) act' (lo) 0/3' ('yo)'Y'Y , 0 XY ,a'/3 , '/'( a l' 1 a 3) ■ (3.4a) 

®XY,al3'y( a 3i a 2> a 3) * (ai , a 2l a 3 ) , (3.4b) 

^AT,a^( a i> ° 2 ; 03 ) — ' * X aa 'Xbb'Y cc 'Ox\f a p^(ai, a 2 , a 3 ) , (3.4c) 

where in eq. (3.4b) charge conjugation is performed first. Additionally we know that each 
operator transforms under a translation in n-direction as 

£ > XY,a0'r( a i + ^,02 + 5a, a 3 + 5a) = exp [i5an- P} C>xK )a/ g 7 (ai, a 2 , a 3 ) exp {-iSan- P] , 

(3.5) 

where P is the momentum operator which acts as a generator of translations. Another 
symmetry of the three-quark operators defined in eq. (3.3) is the invariance under the 
exchange of the quark in the first and the second position or even an invariance under 
exchange of all three quarks in case of the operators containing right-handed or left-handed 
fields exclusively. On top of this the operator is invariant if one simultaneously rescales 
a,; —> A ai and —» n M /A, which we will call scaling property. 


(3.2) 


(3.3) 
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i 

t a ® 

II 

4 *,i 

4 *,i 

dT 

d ? 

1 

il ® t/i 

-1 

-1 

-1 

2 

-1 

2 

1 ® 1 

1 

-1 

1 

1 

0 

3 

( 7 dn ® Tjl 

1 

1 

-1 

2 

-2 

4 

a^ n ® 7 p 

1 

1 

-1 

2 

-1 

5 

® 

1 

1 

1 

1 

0 

6 

ia 9n ® 1 

-1 

1 

1 

1 

-1 

7 


1 

1 

1 

1 

-1 

8 

ct ^ 9 ® 7^ 

1 

1 

-1 

0 

0 

9 

cr^ n ® a^ n 

1 

1 

1 

3 

-2 


Table 2. List of T\ ® ^ 4 = 1 by choice (see comment in the text). We have multiplied 

structures 1 and 6 with a factor of i such that i = rff i for all structures and, thus, r ]^ 0 = 1 . 
In cases where four-vectors are used in the place of Lorentz indices the notation means that the 
corresponding Lorentz index is contracted with the index of the vector; e.g. a dn = a^d^n^. 


3.2 Low energy operators 

Using the previously defined fields ur and ur we can write down the operators, which 
contribute to baryon-to-vacuum matrix elements of three-quark currents at leading one- 
loop level and have correct transformation properties under chiral rotations in the following 
compact form: 


®XY,aP 7(at, 0,2 , <Z3) 


f [dx] E E Fxy( x 1, x 2, *3) 

i j k=1 


r}j,k,XXY 
^S.abc 


(*) 


(3.6) 


where the correct transformation behaviour under translations in n-direction is ensured by 
^ = n M £Xidi and the constraint that + X 2 + X 3 = 1 in 


J~ \dx\- J ° dx\dx 2 dx‘i <5(l - E ®i) 


(3.7) 


where the integrations run from 0 to 1. The JF’s are functions of x\, X 2 , X 3 only and kj is 
given in table 3. The T’s are defined as 


r S Z = (ixT^YC^zF^i^Mn ■ d) d ? , (3.8) 

where T\, V l B , d™ and df can be taken from table 2. The occurring derivatives act on the 
.B’s. We have introduced adequate powers of i$ to have functions J- of mass dimension 2, 
which is compatible with the standard mass dimension of distribution amplitudes. Using 
i$ (which leads to a factor mg in the final result) instead of the baryon mass in the 
chiral limit mo (which would be the standard choice) has the advantage that it allows for 
a straightforward matching of our parametrization to the general decomposition given in 
ref. [ 6 ] and to refs. [7, 8 ] (see also section 4.3). The power of ( n-d ) is chosen such that the 
scaling property is fulfilled. Note that in the chiral-odd sector one can actually write down 
more structures, which have the form or ^ 444 ' H° weve U these structures are 
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j 

jii 

1,5 

B 7 

Bi 

trace 7 

kj 

1 

B S 

t 

t 

1 

3 

2 

B s 

1 

t 

tr {x+} mj 2 

3 

3 

b 5 

1 O 

+ 

>X 

t 

1 

6 


Table 3. In this table we list only terms which contribute to the one-loop calculation of baryon- 
to-vacuum matrix elements of the operator. x+ is defined as x+ ~ tr(x+) /3. This is a convenient 
choice since this combination (in a leading one-loop calculation) vanishes along the symmetric line, 
where m u = rrid = m s . 


not independent. They can be rewritten in terms of using Fierz transformation. 

In order to reduce the T’s to the minimal set given in table 2 one has to use the identity 
cr ,w 75 = ^£ pupa a pa and the fact that it is sufficient to construct structures of positive parity 
(see explanation below eq. (3.15)). Additionally one has to use that multiplying both 
structures Y l A and V l B with a 75 does not lead to a new, independent structure owing to 
the projection with 7 l /r in eq. (3.8). 

The B 's in eq. (3.6) are defined as 


where 


BS,abc ~ ( u x)aa'(uY)bb'(uz)cc'B^ a , b , c , , 

(3.9) 

Bg’abc = (.Byft) aa' (B'i^hb' (-^3 )cc'^a' 6 'c' x trace 7 , 

(3.10a) 

Bg’abc = (B[> )aa/(.B'^ b )bb' {.B^cc'^a'b’c' x trace 7 , 

(3.10b) 

Bs,abc = (B ‘2 )aa' (Bg )bb' (B^ g )cc'£ a 'b'c' x trace 7 , 

(3.10c) 

Bs,abc = ("®2 ^aa’ (B^^bW (^3 )cc'^a'b'c' x trace 7 , 

(3.10d) 

Bg’abc = {B^ ^aa'iB^hh'^B^cc'^a'b'c' x trace 7 , 

(3.10e) 

Bfi’abc = (B'ij)aa'(.B' 2 )bb'(,B 2 \ b )c, c '& a 'b' c ' x trace 7 . 

(3.10f) 


For cases where B 3 2 = B^ we only use B^\ , Bi’ 2 bc and Bl’ :i h and thus kj = 3. The different 
possible combinations of B’s can be taken from table 3. All baryon and meson fields which 
are connected to each other (by a summation over a shared flavor index) have to be at 
the same spacetime position, owing to the fact that the compensator field K is a local 
transformation. However, chiral symmetry actually also allows for the possibility that the 
trace term in B is situated at a different spacetime position as the rest of the operator. We 
consider this possibility in appendix D and show that such a parametrization only differs 
in higher order terms. Note that no structures of the form [B ( s,x+], {B$,x + }, or tr{H, 5 \' + } 
occur in table 3, since they can be reexpressed in terms of the third structure, which 
means that we have only one second order structure (j = 3) that is responsible for SU(3)y 
breaking. Also the operators which describe the behaviour along the SU(3)y symmetric 
line (j = 1 , 2 ) are not linearly independent, but the situation is more complicated in this 
case: since operators of the same class (i.e. same j but different k) are related to each other 





(see eq. (3.30)) one has to take care that the symmetry properties of the operator under 
quark exchange are respected. Therefore, we postpone this discussion to section 3.4. 

There are no covariant derivatives acting on the baryon field within the B : s. In ap¬ 
pendix D we show that they can always be traded for derivatives acting on the whole 
structure plus higher order contributions, which can be neglected. This fact will turn out 
to be very convenient for calculating loop contributions, since the derivatives acting on the 
complete structure do not lead to additional loop momenta in the integrals. 

The effective operator given in eq. (3.6) already transforms correctly under chiral ro¬ 
tations and translations along the light-cone vector n. It also fulfills the scaling property. 
The remaining symmetry properties given in section 3.1 will now be implemented by con¬ 
straining the functions T. We consider 


and 


R j,k,XYZ _p_ 
<5, abc 




nj,k,XYZ d t V n jjj,k,XYZ 
B S,abc —* ~ C SS'B J 5 , jabc , 


(3.11a) 

(3.11b) 


r i,XYZ P 
afY/S 

v i,XYZ c t f 
afY/d 


'Vr,i ('Yd)aa' ('Yo)/30 , ('7o)'y'y' ^-'a'j3'YS' (To)<5'<5 


-Vr P i C CU'C^CV rcw 


,i,XYZ 


(3.12a) 

(3.12b) 


Eqs. (3.11b) and (3.12b) yield (together with eqs. (3.4b) and (3.6) and since r^ PC = 1) 


vxy ) r : 


XY ’ 


(3.13) 


which would mean that the J-’s are real-valued. However this argument relies on the 
assumption that one gets no additional phases from charge conjugation of quarks and 
baryons, which is not necessarily true. If we allow for such additional phases the above 
equation has to be generalized to 

{^e i9 Y =^ k e ie , (3.14) 


where we have an additional overall phase which is equal for all distribution amplitudes. 
However, this additional phase is unphysical and can be dropped. Eqs. (3.11a) and (3.12a) 
yield (together with eqs. (3.4a) and (3.6) and since r/f = 1) 


sT~'i . 7 .k st~t . 7 ". k 

B xy = ~ T XY 


(3.15) 


7 i,j,k 


Therefore we only have to differentiate between chiral-even - —j LL 


tri.j.k _ T-i,j,k i 

J~ T T — J~ ^ven Q'licl 


chiral-odd 


Notice that we have chosen to only construct structures 
T^ ® T# which have positive parity. The negative parity structures, which one can obtain 
by multiplying all T# with a 75, would lead to the same operators since eq. (3.15) then 


would yield an extra minus sign. 
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3.3 Symmetry under exchange of quark fields 

In this section we use the symmetry of the original three-quark operators under exchange 
of quark fields with the same handedness to reduce the number of amplitudes. Using the 
constraint that the operators have to be equal under exchange of the first and the second 
quark yields 


3 = 1,2: 

f&Hx 1 , 

x 2 , 

,x 3 ) = 

~ r lr A ,1 -^xy ( x 2,x\,x 3 ) , 

(3.16a) 


1 , 

x 2 , 

,x 3 ) = 

-VrA'iFxY( x 2’ x i’ x s) , 

(3.16b) 

3= 3: 

T i,S,1 (x 1 

x 2 , 

,X 3 ) = 

_r ?r A ,i“^vy 4 ( ;E2 ’ x i’ x 3) , 

(3.16c) 


** XY 

x 2 , 

,X 3 ) = 


(3.16d) 


'T’i, 3,3 / 
•'XY 

x 2 , 

,X 3 ) = 

-r/p .4 Yxy{ x 2 -, x \, x 3 ) • 

(3.16e) 


In the chiral-odd sector one now uses these relations to eliminate (if 3 = 1,2) and 

/ 5 / fi . Additionally we can use that 

(T^r AlxC) 1 p{^ x ^ Es)a5 = 0 , if X ±Y and T A e {1, 75 , 07 ^} . (3.17) 


Using Fierz transformation this leads to 

r 3,XAT _ r 4 ,XXY lp9 ,XXY 

aft'yS afl'yS 2 ’ 

r 5 ,XXY _ n 
1 a/3^/6 ~ U > 

p 6 ,XXY _ r 4 ,XXY _ r 7 ,XXY 
afl'yS a/3'yd afi'yS ’ 

if A + Y. Therefore we have the freedom to choose 


(3.18a) 

(3.18b) 

(3.18c) 


T-3j',fc 

“'odd 


(x 1 ,X2 : X 3 )=J*dd 


(xi,X 2 ,X 3 ) = 


•6 ,j,k 
odd 


(xi,x 2 ,x 3 ) = 0 . 


(3.19) 


In the chiral-even sector the projection with Jl/r leads to similar constraints. The coun¬ 
terpart of eq. (3.17) reads 


{lx?AlxC)~tp(ixTB)cx6 = 0 , if T a e { 7 ^ 7 ^ 75 } ■ 

With a Fierz transformation one obtains 

r 7 ,xxx _ _ r 4,XXX 1 r 5,XXX r G,XXX 
a/3^5 2 cx(3^8 ’ 

r 8,AA.V _ _]_ t 5.XXX 
a/37<5 ^ a /375 ’ 

r 9,XXX _ n 
1 a/3 7 5 ~ U ' 

Therefore, we can choose 

Flven (x 1 ,X 2 ,X 3 ) = Tl^{xi,x 2 ,x 3 ) = F%$£(xi,x 2 ,x 3 ) = 0 . 


(3.20) 

(3.21a) 

(3.21b) 

(3.21c) 

(3.22) 
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The operators containing left-/right-handed quarks exclusively also have to be invariant 
under an exchange of the first and the last quark. Performing a Fierz transformation and 
using the identities given above we find 



r i,xxx _ 

r i', XXX i> 
2^ cx/3 7(5 ^ 

'i 




i-l 




The matrix c is given by 


/ i 

2 

0 

1 

2 

3 

2 

0 

2 


0 

1 

2 

0 

0 

6 

1 

2 


0 

0 

1 

0 

0 

0 

c = 

1 

2 

0 

1 

2 

1 

2 

0 

1 

2 


0 

1 

8 

0 

0 

1 

2 

1 

8 


l 0 

0 

0 

0 

0 

1 


(3.23) 


(3.24) 


By the use of this relation the symmetry property of the operator under exchange of the 
first and the last quark translates to the following constraints on the amplitudes: 


3 = 1,2: 


3 = 3: 


^ XX 

(xi,x 2 ,x 3 ) = 

6 

- E ■=“ 

j'=l 


»2, 

x l) , 

(3.25a) 

‘'XX 

(xi,x 2 ,x 3 ) = 

6 

- E 

j'=l 


»2, 

x l) , 

(3.25b) 

3 

^ XX 

(xi,x 2 ,x 3 ) = 

6 

- E 

j'=l 


»2, 

x l) , 

(3.25c) 

3,1 

(xi,x 2 ,x 3 ) = 

6 

- E 

j'=l 

TYX Y x 3: 

( X 2: 

, x l) , 

(3.25d) 

T-i,3,2 
TYX 

(xi,x 2 ,x 3 ) = 

6 

- E ■=“ 

j'=l 

TYX V x 3- 

i®2: 

, x l) , 

(3.25e) 

/j^i, 3,3 
‘'XX 

(xi,x 2 ,x 3 ) = 

6 

- E c“ 

j'=l 

~ r XX Y x 3- 

1^2: 

,ar) , 

(3.25f) 

3,4 

•'XX 

(xi,x 2 ,x 3 ) = 

6 

- E # 

i'= 1 

TYX V x 3- 

J 2: 

,®l), 

(3.25g) 

^*,3,5 

•'xx 

(xi,x 2 ,x 3 ) = 

6 

- E 

i'= 1 

V'’ 3 ’ 3 ^, 

TYX V x 3- 

1^2: 

, x l) , 

(3.25h) 

,3,6 

•'xx 

(xi,x 2 ,x 3 ) = 

6 

- E 

TYX V x 3- 

,X 2: 

,Xl) . 

(3.25i) 


i '=1 


Using these equations one finds for the operator with j = 3 that one can eliminate all 
amplitudes apart from J'^x , by using the following relations recursively: 


?xx ( x i, x 2 , x z) = VrA,i E ° U , 


t ii r st,3,3 / 


(3.26a) 


i '=1 
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/T72,3,3 

^XX 


(XI,X 2 ,X 3 ) 


11° 
r lr A a 


6 

E c " 


r 


,3,1 


XX 


i'= 1 


(x 3 ,x 1 ,x 2 ) , 


^\ 4 (x 1 ,x 2 ,x 3 ) = -VrA^x 3 x(x 2 ,x 1 ,x 3 ) , 

^(11,12,2:3) = -VrA^xx ( x 2 ,xi,x 3 ) , 

Fx 3 x( x i, x 2, x 3) = -VrA ti F l xx( x 2,xi,x 3 ) , 


For the operators with 7 = 1,2 we can eliminate 

E\4 2 ( x i’ x 2,Z3) = -VrA ti ^xx( x 2,xi,x 3 ) , 

^X^bl'2^3) = - E C M J’^ 1 (X 3 ,X2,X 1 ) , 
i'=l 


(3.26b) 

(3.26c) 

(3.26d) 

(3.26e) 

(3.27a) 

(3.27b) 


and additionally 


JF 


x J x ( x !, x 2, x 3 ) = ^^(xi, x 2 , x 3 ) + Fxx(xi, x 2 , x 3 ) - 2 ^y(xi, x 3 , x 2 ) 


+ ^^ 1 (xi,x 2 ,x 3 ) , 


^( Xl ,x 2 ,x 3 ) 

^^. 1 (xi,x 2 ,x 3 ) 

J r ^^ 1 (xi,x 2 ,x 3 ) 


= -4J r ^ 1 (xi,x 2 ,x 3 ) + 8^ 7 ^ 1 (xi,x 3 ,x 2 ) + JF^ 1 (xi,x 2 ,x 3 ) , 
= -^^ 1 (xi,x 3 ,x 2 ) , 

= -F*xx( x 1, x 3, x 2) ■ 


(3.28a) 

(3.28b) 

(3.28c) 

(3.28d) 


From the fact that the local operator at the origin, where a\ = a 2 = a 3 = 0, is indepen¬ 
dent of the light-cone vector n one can deduce constraints for the zeroth moments of the 
distribution amplitudes 

J[dx]J^Y{x i,x 2 ,x 3 ) = 0, for i = 1,3,4,6,7,9 . (3.29) 

3.4 Elimination of linearly dependent structures 


To avoid overparametrization we will now annihilate linearly dependent structures of those 
given in table 3. Considering all possible three-quark operators and all baryons from the 
octet, one finds (for j = 1,2) that only two out of the three structures B J s ’^ bc , B J s ’ 2 bc and 
Bfi’^bc are linearly independent, since one has 


0 “ B i\abc + 


Rh 2 + R.?’ 3 


(3.30) 


In the chiral-odd sector we can use this relation to replace B^ ab = ~Bf l fibr - B x \ , which 
is equivalent to the replacement 


^ 1 (*l,®2,* 3 ) 

J^(X U X 2 ,X 3 ) 

?rid( x U x 2, x 3) 


^odd (®i,®2,a: 3 )=^d 1 
Kdd( x ^ x 2, x 3) = Kdd 
Kdd( x ^ x 2’ x 3) = 0 • 


(xi,x 2 ,x 3 )-^*^ 3 (xi,x 2 ,x 3 ) , 
(xi, X 2 , x 3 ) - JF*’^ 3 ( Xl ,x 2 , x 3 ) , 


(3.31a) 

(3.31b) 

(3.31c) 


Using eqs. (3.16a) and (3.16b) one finds that the new functions have the same symmetry 
properties as the old ones. Therefore we can choose 

3\,m( Xi ’ X2 ’ x z') = 0 > J = 1 ’ 2 ’ ( 3 - 32 ) 
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in accordance with symmetry properties and without loss of generality. In the chiral- 
even sector the situation is different since the amplitudes are already constrained by the 
symmetry under exchange of the first and the third quark. An elimination of one structure 
in favor of the two others would therefore not lead to a simplification. Instead one just 
obtains a reparametrization of the problem for which it would be hard to implement the 
symmetry properties under exchange of the first and the last quark. 

4 Calculation at leading one-loop order 

In this section we describe the leading one-loop calculation. In section 4.3 we explain how 
we have matched to the standard DAs defined in ref. [6]. 

4.1 Meson masses and the Z-factor 

We work in the limit of exact isospin symmetry, where m u = m^ = mi. Using the standard 
leading order meson Lagrangian (see e.g. [16, 21]) one finds for the meson masses the 
standard Gell-Mann-Oakes-Renner relations 


ml = 2 B 0 mi = m- =li2j3 = 2 B 0 (m q - 8 m t ) , 

(4.1a) 

m 2 K = Bo(m t + m s ) = mf =4 7 = B 0 (2m q + 5mi) , 

(4.1b) 

m^ = —! ~(2mi + 4m s ) = m* =8 = 2 B 0 (m q + 5m t ) , 

O 

(4.1c) 

1 

m q = -(2 mi + m s ) , 

(4.2a) 

5mi = fh q - mi , 

(4.2b) 


and Bq is the LEG proportional to the quark condensate in the chiral limit. As additional 
ingredient we need the first order meson-baryon Lagrangian, which we take from ref. [22] 
(this version differs from refs. [18, 21] by a minus sign in the terms containing D and F in 
order to be consistent with the standard sign convention “ D + F > 0): 

4/b =tv{B'y ti iD^ l B}-m 0 tr{BB} + ^ tr {B'y li 'y 5 {u fl , £}} + ^ tr {Bj^ 5 [u M , B]} . 

(4.3) 

For our calculation we need the baryon-meson-baryon vertex for an incoming baryon B , 
an outgoing baryon B 1 and an incoming meson (the fc-th one in the Cartesian basis) with 
momentum q , which is given by 

-±^ 5 tr {K T B ,(D{\ k , k b } + F[\ k , k b ])} . (4.4) 

2c 0 

The self-energy to third chiral order is given by the sum of the irreducible diagrams shown 
in figure 1 (where external legs are to be amputated) multiplied with an i. The contribution 
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q 


<■> 

(a) 


p Ax p-q Ax p 


<■> 


0 


(b) 


Figure 1. Feynman diagrams needed for the calculation of the self-energy. 


of diagram (b) (in figure 1), which is relevant for the calculation of the Z-factor is given by 
i* (b) = 3g B , 7T f ( 171 ^, 1710 ,tf) + 4g BtK f(m K ,m 0 ,y>) + g B)ri f(m v ,m 0 ,]/>) , (4.5) 

where 

r(i)/ 


-1 


f(m,mo,$) = -jjA(p 2 - (m,rn 0 ,p) + (f + mo)(l 01 (mo,p) - m 2 In(m,m 0 ,f))) . 

(4.6) 

The loop functions I^i and 1^ are defined as in ref. [23] and the coefficients are given by 


(D + F) 2 , 
4. 


9n, 

5s, 

V 

g~,n = (D-F) 2 , 


~(D 2 + 6F 2 ) 
9 V ’ 


5 A, t 


4 9 
-D 2 
3 


g N ,K = \d 2 - DF + ^F 2 
0 2 

5s,x = D 2 + F 2 , 

9 e,k = ^D 2 + DF + ^F 2 
0 2 

9 a,k = \{D 2 + 9F 2 ) , 


These constants fulfill the constraints that the sums 


3 

4. 


9N,rj = ~(D - 3 F) 2 

? 

4 n 2 

5s, v = £ D > 
gz, v = \{D + 3F) 2 

1 

4 2 

9A,rj “ 2 ^ * 

(4.7) 

9 F 2 ) , 

(4.8a) 

9 F 2 ) , 

(4.8b) 


are independent of the baryon/meson species. This yields similar baryon masses and Z- 
factors along the line of equal quark masses and is a consequence of SU(3)y symmetry. For 
a detailed study of baryon masses under symmetry breaking see [24] . The square root of 
the Z-factor needed in our calculation is given by 


\[~Zb - 1 + -£b > 


(4.9) 


where the prime indicates taking a derivative with respect to f and substituting ]/> -*■ m B , 
while 

Y! B {m q ,8mi) = 3g Bt7T f'(m n ,m 0 ,m B ) + 4 g B , K f{mK,mo,m B ) + g B ^f{m r) ,mQ,m B ) 

= T/*{fh q ) + AE' B (m q ,6mi) , 

(4.10) 
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with 


£'*(m g ) = Z,' B (m q ,0 ) = (3 g B ^ + ^9b,k + 5 s^)/'(m^,m 0 ,ms) 

= ^(5T > 2 + 9F 2 )f\m^,m 0 ,m B ) = ^(5T > 2 + 9F 2 )f\m* n ,m* b ,m* b ) 


(4.11a) 


A'Z' B (m q ,5mi) = Y 1 ' B (m q ,5mi) - F>' B (m q , 0) 


qi^uvi) - ^gynvq,vuvL j z-igyuuq, 

= 3 gB , 7 tf '( rn n , m B , m B ) + ^9b,k Y {rn K ,rn. 0 ,m. B ) + gB^f (rn v ,m, 0 ,m, B ) 
~ \^ D2 + 9i ?2 )/'(m^,m 0 ,m s ) 

= ^9B,-Kf'(rn n ,m^,ml) + Ag BiK f'(mK,ml,ml) + g BtV f'(m v ,ml,ml) 


)(5D 2 + 9F 2 )f'(m* m ,ml,m* b ) , 


(4.11b) 


where = m* m j h (m q ) is the meson/baryon mass along the symmetric line ( 5mi = 0 ). 

The dotted equal sign = means equal up to terms which are of higher order than our 
level of accuracy (which is second order in chiral power counting). For explicit results see 
appendix A. 


4.2 Baryon-to-vacuum matrix elements of three-quark operators 

In this section we describe the actual loop calculation. From a simple power counting 
argument one finds that at leading one-loop level the only contributing graphs are the ones 
shown in figure 2 . One easily observes that the second order operator insertions only occur 
without additional mesons. Therefore we only have to compute the vertices where a single 
baryon couples to the operator. Contributions with additional mesons only occur for the 
leading order operator insertion (j = 1). For the BChPT calculation mainly the structure 
Bfabc * Z re l evan t- Graph (d) of figure 2 is an exception because the extra 75 from 
the baryon-meson-baryon vertex has to be canceled with a 75 from the Dirac structure 
of the operator. The calculation gets simplified considerably if one uses the fact that (by 
construction) the z with k 7 1 can be obtained from the case k = 1 by a permutation 

of indices: 


R j,2,XYZ _ R j,i,vzx 

^S^abc ^Sibca 


B^XYZ = B 3,kZ XY 
d,abc o,cab 


R jA,XYZ _ _ Rj: l.YXZ 
^d^abc ^Sibac 


R j^XYZ _ _ R jXXZY 
^d^abc ^S^acb 


r>jfi,XYZ _ _ n jXZYX 
^d^abc ^S^cba 


(4.12) 


which means that we actually only have to calculate the case k = 1. Defining 


f+1, for X = 
j-1 , for X = L 


we can write down the relevant operator insertions in a quite economic way: 


(4.13) 


j^l,l,XYZ 

S,abc 


(*) 


= {kbY 


_a'bc -ip • 2 : 


Sse 


Be( P ) 


(4.14a) 
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p 


(a) 



P 


(b) 


q 



q 


/ 


\ 


/ 


\ 

\ 


\ 


/ 

/ 

p 

\ 

0 





(c) 


p 


(d) 


Figure 2. Feynman diagrams needed for the calculation of the baryon-to-vacuum matrix ele¬ 
ments. The squares depict the operator insertions given in eqs. (4.14), the circle stands for the 
vertex from the meson-baryon Lagrangian given in eq. (4.4) and the dashed/solid lines represent 
mesons/baryons. Diagram (a) has to be multiplied with \fZ. However, one knows that at higher 
orders all of the diagrams will receive a \fZ contribution, which can be used as an argument in 
favor of the factorized version of our results (see eq. (5.18) in section 5). 


tj1,1,XYZ 

8,abc 


(*) 


Be{p~q)<t> k {q) 


(-l)x(A k K B ) aa '5 w 5 cc ' + (-l) Y (n B ) aa \X k ) w X™' (4.14b) 

Zrr\ L 


+ (-l)z(KB) aa 'S bb '(\ k ) CC ' 


a'b'c' -ip - z 


Sse , 


B 


bb xyz (A 

8,abc W 

-1 


16 F$ 


Be (p-qi -|? 2 ) <t> k (qi) <t> 1 (<? 2 ) 

[({A fc ,A l }K B ) aa '5 W 5 CC ' + (n B ) aa \{X k ,X l }) bb 'S cc ' + ( K B ) aa '5 bb \{X k ,X l }) cc ' 
+ 2(-l) x (-l) Y (X k K B ) aa '(X l ) bb 'S cc ' + 2(-l)x(-l)z(X k K B ) aa ’d bb ’(X l ) cc ' 

+ 2{-l) Y {-l)z{K B ) aa \X k ) bb \X Y + 2{-l) x {-l) Y {\ l K B ) aa '{\ k ) bb '5 cc ' 
+ 2{-l) x {-l)z{X l K B ) aa ' 5 bb \X k ) cc ' + 2{-l) Y {-l) z (K B ) aa \\ l ) bb \X k ) ccr 


w jalb'd —ip ■ z r 
x 8 e y OSe • 


The second order tree-level operator insertions read 


d 2,1 ,xyz ( n 

B s’abc ( 2 ) 


= 450^0^ tr {A4} (n B ) aa '£ a ' bc e~ ip z 6g e , 


B 


3,1,XYZ 
S.abc 


(*) 


S«(p) 


B«(P) 


45 0 mo 2 (K B ) aa, 7W 6b, e a ' b ' c e- ip ‘, 


(4.14c) 

(4.14d) 

(4.14e) 
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where JA = JA -tr {JA}/3. After performing the loop calculation one finds that the results 
can be expressed as 


(0| £ > RR,aP'y( a h a 2’ a 3) + ^LL■> a 2 i a 3)\B(p, s)) - 

- f WC’, 

(0\O l R[ ita p^( a i,a2, a 3) + Ofy ^(( 11 ,( 12 , a 3 )\B(p, s)) = 

- f [*]C"' £rt “ErS«?(w)4*( I1 ,» ! ,B), 

J i 

where uf(p,s) is the baryon spinor, 

pi,even _ j^i,RRR -pi,LLL 
aP'yS a(3^5 Oifi^yS ’ 

■pi,odd y\i,LLR "pi,i?i?L 
oifi'yS a(3~/5 afljS ’ 


(4.15a) 


(4.15b) 


(4.16a) 

(4.16b) 


and 


hstve n (® 1, ®2, ®3) = E ( "R,RRTlRven 0*4, , X\ 3 ) , 


j,k 


li.abc / \ j.k.abc rri.j.k/ \ 

h B,odd( X ^ X 2, x 3) = E^LLR^odd Ol^2,Z 3 ) • 

j,k 


(4.17a) 

(4.17b) 


The coefficients c> B xyz inherit the property that the ones with k ± 1 can be obtained from 
the case k = 1 by a permutation of indices: 


i,2 ,abc _ i,l,6ca 
C B,XYZ ~ C B,YZX ’ 
i, 5 ,a 6 c _ i,l,acb 

C B,XYZ ~ ~ C B,XZY ’ 

For those with fc = 1 we find 


i, 3 ,a 6 c _ i,l,ca6 

C B,XYZ ~ C B,ZXY J 


i, 4 ,a 6 c _ i,l,bac 

C R.XYZ ~ ~ C B,YXZ ’ 


i,6,a6c 

C B,XYZ 


i,l ,cba 
~ C B,ZYX * 


1 ,1, abc _ 1,1, abc 
-B,XYZ ~ C B,XYZ 


+ C 


(a) 


1 ,1,abc 
B,XYZ 


+ C 


(c) 


1 ,1,abc 
B,XYZ 


(d) 


1 ,1, abc 
"B.XY Z 


o l,l,a6c 

2 B,XYZ 


= Yz^(K B ) aa 'e a ' bc , 


(4.18) 

(4.19a) 

(4.19b) 


(a) 


-1 


1 ,1, abc 
C B,XYZ 


~ 8F o T[( xkxk KB) aa 's bb 's cc ' + {B R ) aa \\ k \ k ) hh '8 CC ' + (K B ) aa '8 bb '(X k X k ) cc ' 

(c) 0 +2(-l) x (-l)y(A fc «s) aa '(A fc ) bf,, 5 cc ' 

+ 2(-l) x (-l) z (A fc « s ) aa '^'(A fe ) cc ' 

+ 2(-l) y (-l) z (ic B )“ a '(A fc ) 66 '(A fc ) cc ']e“' 6 ' c 7 1 o(mO , 

(4.19c) 

x (/io(m fc ) + (ra| -rUfl)/ii(m fc ,m 0 ,m B ) -m B (m B + mo)/ 1 ( } ) (m fe ,mo,m s )) . 

(4.19d) 
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In the contribution from graph (d) commuting 75 from the vertex with the Dirac structure 
in the operator yields ^(- 1 )^* = -1 (compare table 2 ). In operators of type Oxr the 
75 has no effect owing to 7^75 = 7 r. The relative sign in the vertex in operators of type 
Oxl i s compensated by 7^75 = - 7 l- Therefore the result only contains structures of the 
form given in eq. (4.16). This is no coincidence but has to happen in order to obtain a 
result that behaves correctly under parity transformation. For the second order tree-level 
contributions we find 


o 2,l,a6c 

~B,XYZ 


aa' _ a'be 


4 B 0 m 0 2 tr {M} (kb)“ s' 


(4.19e) 

4 )xyz = 4B 0 m- 0 2 (K B y a 'M W e a ' b ' c . (4.19f) 

Using eq. (3.2) the matrix element of the complete three-quark operator reads 

{0\q“(ain)q b p(a 2 n)q°(a 3 n)\B(p, s)) = 

[dx]^ " ' PEfcIfc “ fc Y,{ T aJj 5 h B%n(xi,X 2 ,X 3 ) + T^h^^Xu X 2 , X 3 ) 
i 

+ T yap8 h B, odd( X 3, X 1 ,X 2 ) + T^ aS h b,odd^ Xl > j U <5 (P» S > ' 

(4.20) 


/> 


4.3 Projection onto standard DAs 

In this section we relate our parametrization of the baryon-to-vacuum matrix element, 
which was guided by the behaviour under chiral rotations, to the general decomposition 
given in ref. [ 6 ], which is more convenient for daily use. To do so we have contracted 
both our result (eq. (4.20)) and formula (2.3) of ref. [ 6 ] with Dirac structures of the form 
0 U® r It is sufficient to use structures where Lorentz indices are either contracted 
between T ' 4 and T B or with the light-cone vector n or the momentum p. Afterwards 
we have used the identity fu B {p) = m B u B (p ) and have matched the prefactors of the 
remaining Dirac structures (75 and ^ 75 ). Using twist-projection, we obtain the results for 
the distribution amplitudes S B , P B , Af, V B and T B which are independent of the scalar 
product n-p, due to the scaling property described in section 3.1 (For the details of the 
twist-projection we refer to [ 6 ]). We have collected these lengthy matching relations in 
appendix E. The amplitudes have the following symmetry properties under exchange of 
the first and the second variable 


where we use 


Sf(xi,x 2 ,x 3 ) = -(-1 ) b S b (x 2 ,x 1 ,x 3 ) , 
P b (x 1 ,x 2 ,x 3 ) = -{-1) b P b (x 2 ,x\,x 3 ) , 
Af{xi,x 2 ,x 3 ) = -(-l)sAf (x 2 ,xi,x 3 ) , 
V b (x u x 2 ,x 3 ) = +(-l) B V B (x 2 ,x 1 ,x 3 ) , 
T b (xi,x 2 ,x 3 ) = +(-1) b T b (x 2 ,xi,x 3 ) , 


f+1 , for B ± 
[-1 , for B = A 


(4.21) 


(4.22) 
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for brevity. To obtain these nice symmetry properties one has to choose the flavor content 
in the operator as p = uud , n = ddu, S + = uus , S° = uds , = dds , H° = ssrt, = ssd, 

A = uds, where the order of the flavors is relevant. The different sign for the A originates 
from the antisymmetry of the isospin singlet state. 

5 Results 

In this section we present our results and provide a definition for DAs that do not mix under 
chiral extrapolation. In section 5.3 we work out an explicit parametrization of baryon octet 
DAs, where we follow the approach presented in refs. [9-11]. 

5.1 General strategy and choice of distribution amplitudes 

We will split up every distribution amplitude in the following way: 

DA (fh q , Smi ) = DA(m q , 0) + (DA(m,j, Smi ) - DA(m ? , 0)) 
e DA*(?h g ) + ADA(fh q ,5mi) , 

DA*(m 9 ) = DA*(0) + (DA *(m q ) - DA*(0)) 

EDA° + ADA*(m 9 ) , 

where the main idea is to use the second formula to parametrize everything in terms of 
the DAs at the symmetric point, which are measurable on the lattice as opposed to the 
amplitudes in the chiral limit. Lattice simulations where the mean quark mass is fixed at 
its physical value while 6mi is varied are already available for hadron masses and some 
form factors [25-27]. Corresponding simulations for the baryon octet DAs treated in this 
work are in progress. This strategy has the additional advantage that one gets rid of the 
parameters that describe the behaviour under variation of the mean quark mass. For the 
presentation of the results it turns out to be convenient to write down the second order 
tree-level and the loop contribution separately. We define for all baryons 

ADA = ADA loop + 5m ADA tree , (5.2) 


(5.1a) 

(5.1b) 


where 


5m = 


4Bo5mi 



(5.3) 


Then we use the fact that we can rewrite ADA in terms of m£ and DA* using the cor¬ 
responding expansions in fh q . For a specific set of DAs, which do not mix under chiral 
extrapolation (see below), this allows us to rewrite the loop contribution as the DA along 
the symmetric line multiplied with a loop function / such that the results have the form 


DA(m ? , 5mi) = DA*(m ? )(l + /) + 5m ADA tree . (5.4) 
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By virtue of SU(3)j symmetry we find the following relations between DAs along the line 
of symmetric quark masses m u = rrid = rn s : 


2T 1/ e(x l ,X 2 ,X 3 ) = ( -1 ) B [ - A fy* ] (x 1 , X '3 , X' 2 ) 

+ [Vye -^?/ 6 ](* 2 ,® 3 ,® i ) , 

[Till + T,% + - Pi / ;\(x^x 2 ,x 3 ) = [F 2 %* - A^\(x 2 ,x 3 , Xl ) 

+ [Vyt -Afil](x 3 , X l,X 2 ) , 

2 T 2 / 3 (xi,X 2 ,X 3 ) = [T 3/4 - T 7 y g + .S) j 2 + P 1 y 2 ](^3)^1 1 ^ 2 ) 

+ [^4*-^8*+'5f/; + ^ 2 1(^3,X 2 ,X 1 ) . 


(5.5a) 


(5.5b) 


(5.5c) 


Note that we do not impose these relations. They are automatically fulfilled by our calcula¬ 
tion (loop contributions included). For the nucleons these relations are fulfilled exactly also 
for 5mi 4- 0 owing to isospin symmetry (again this is also true for the loop contributions), 
which was already shown in ref. [ 6 ]. If we were only interested in the SU(3) j symmetric case 
(or in nucleons only), it would therefore be enough to define the independent amplitudes 
as 


$f / 6 (xi,s 2 ,® 3 ) = \Vy 6 - A^ /6 \(x 1 ,x 2 ,x 3 ) , (5.6a) 

$f / 5 (xi,x 2 ,a; 3 ) = [Vq 5 - ^ 5 ](® 1 ,® 2 ,® 3 ) , (5.6b) 

^ 4 / 5 ( x u x 2 ,x 3 ) = [V 3 / 4 - Ay 4 ](x 1 ,x 2 ,x 3 ) , (5.6c) 

“ 4 / 5(^15 x 2 i x 3 ) = [T 3/4 — T 7 /g + 5 i / 2 + Pi/ 2 ](a;i, x 2 , x 3 ) , (5.6d) 


where the and 'S’f describe the coupling to chiral-odd operators, while the Ef describe 
the chiral-even sector. The subscript indicates the twist. As it turns out the amplitudes 
<&f, 'hf and Ef are not yet the optimal choice for a description of the complete baryon 
octet, since they mix under chiral extrapolation. Additionally one finds that it is very 
convenient to use differing definitions for the A, which we choose in such a way that the 
DAs of the A coincide with the DAs of the other octet baryons in the limit of equal quark 
masses. Therefore we define 

^±, 3 l§{ x l-> x 2-, x 3) = “([^ 1/6 “ 

^±, 4 / 5 (* c i > X2 ’ x s) = c s([^2/5 ~ ^2/5] 

E f,4/5( X '1, X '2,T3) - 3(-l) S c|( [Tf /4 -l 

±[T§i- 

where 

+ J1 , for B * A 

B l\/i. for B = A ’ 

Being interested in SU(3)j violation one can not use the constraints given in eq. (5.5) and 
therefore one needs six additional DAs. Our choice are (up to differing prefactors for the 


,x 2 ,x 3 ) A 4/6 \(x 3 ,x 2 ,xi)) , (5.7a) 

,X 2 ,X 3 ) ± (-l)s[v r 3 y 4 - Af /4 ](x 2 ,x 3 ,xi)) , (5.7b) 
7/8 + 5 f/2 + p \%]( x i > x ‘2,x:i) (5.7c) 

- T i/8 + S i/ 2 + Py 2 ](xi,x 3 ,x 2 )) , 


C B 


11 , for B * A 
I -y/E , for B - A 


(5.8) 
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A and exchange of variables) the left-hand sides in eq. (5.5) since they coincide with the 
DAs in eq. (5.7) in the SU(3)j symmetric limit. We define 

n f /6 (xi,x 2 ,a; 3 ) = c B (-l) B Ty 6 (xi,x 3 ,x 2 ) , (5.9a) 

n 4 / 5 ( Xl ,x 2 ,x 3 ) = c b [T 3/4 + T 7/8 + 5i /2 - P 1 / 2 ](x 3 ,x lt x a ) , (5.9b) 

^A/ 5 ( X UX 2 ,X 3 ) = 6c B Ty 5 (x 3 ,x 2 ,xi) , (5.9c) 

where the Ilj describe the chiral-odd sector, while the T, describe the chiral-even part. For 
each octet baryon the standard DAs can be decomposed into the amplitudes defined in 
eqs. (5.7) and (5.9) (see appendix B). We find that the DAs for different nucleons, S’s and 
S’s are related to each other exactly by isospin symmetry. Therefore we define 

DA^ e DA P = -DA n , (5.10a) 

DA S e DA S = -DA S+ = \/2DA s0 , (5.10b) 

DA S e DA h ° = -DA H ~ , (5.10c) 

and give the results only for DA A DA S , DA = and DA a . In the SU(3)j symmetric limit 
all these DAs can be related to those of the nucleon: 


> N * = (j> s * = = <j> A * = n^* = n s * = n s * 

+ ,2 +,^ + ,2 + ,2 2 l l 5 

(5.11a) 

< N * = <j) s * = = $ A * = n A * 

— ,2 — ,2 — ,2 2 J 

(5.11b) 

'TV* ^ A* ^ N * 'Y'S* 

'+,2 — ‘—'+,2 — l_ '+,2 — 1—1 +,2 — -*-2 ~ ^ i ~ ^ i t 

(5.11c) 

>7V* _ 77“* 77-A-* 'V'A* 

‘ — A — '—'—2 — 1 —'—,2 — 1 —'—A — 2 

(5.lid) 


5.2 Minimal parametrization of baryon octet distribution amplitudes 


The choice of DAs presented in the previous section allows us to write down our results in 
a very compact form: 


^(l + iAE^AfffJ + imA®^, 

(5.12a) 

H^l + lAEk + AgfJ+am AS® , 

(5.12b) 

<&l fl)i (l + ±AEk + Asg)+<Sm AUf , 

(5.12c) 

S ±s,*(l + |ASb + Agf ) +5m ATf , 

(5.12d) 


where “±b” stands for “+” if B ± A and for if B = A. The second term in the equations 
above originates from quark mass insertions, while the first term (or, to be more precise, 
AT,' B and Ag® A ) is generated by meson loops and contains chiral logarithms. Owing to 
our choice of DAs the functions Ag^A, which are listed in appendix A together with ATi' b , 
do not depend on the twist of the amplitude. Ag BA and ATi' b vanish for equal quark 
masses {5m = 0). The nontrivial dependence on the mean quark mass of the distribution 
amplitudes i an d S* - is presented in section 5.4. The amplitudes describing the tree- 
level contribution to the SU(3 )f symmetry breaking are not completely free. It holds for 
all distribution amplitudes 


ADA 5 = -ABA N - ADA S . 


(5.13) 


Furthermore, the amplitudes Allf and AYf can be expressed in terms of A i and AE^: 


Anf = a<&* , 

Anf = ~A^ )i -^A^ ii , 
Anf = -^A 4 >(V - ^Ad>f, , 


ATf = AE^, (5.14a) 

ATf = -iAE^-^AS^, (5.14b) 

AYf = -^ASf, : - ^AHfj , (5.14c) 


which means that the Ilf and Y f are completely fixed by the other amplitudes. The 
divergencies of leading one-loop order contained in AE^ and A< 7 jf A can be canceled by the 
introduction of counterterms 


A<h 
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±,i 


m, 


:2 B 
• L $± 


24 F? 


K,iL + A$. ? 


B, ren. 


(/r) , AS 


m, 


±,i 


2 B 

%E1,L + AS B f n '(g) 


24 F 2 




±,l 


(5.15) 


where F contains the divergence and the typical constants of the modified minimal sub¬ 
traction scheme (see eq. (A.5)). F* is the meson decay constant in the SU(3)j symmetric 
limit. The coefficients cf A are given by 

c$± = -9 (F> 2 + 10 DF - 3F 2 ) - 23 t 24 , cf = -9(F) 2 + 10 DF - 3F 2 ) + 9 , 
c| ± = 18(F) 2 - 3F 2 ) + 10 ± 12 , c| = -cf = 18(F) 2 - 3 F 2 ) - 18 , 

4± = -18(F) 2 -3F 2 )-26±12 . (5.16) 


Note that we give no values for c^ ± and eg, since the renormalization of the corresponding 
amplitudes is already fixed via eq. (5.13). The renormalized amplitudes acquire a depen¬ 
dence on the chiral renormalization scale fi, which exactly cancels the scale dependence of 
the leading chiral logarithms: 
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B, ren. 


(m) 


-1 m, c; 


2 R 


•$± 


(4vr ) 2 24F 2 




d 

on 


<B, ren. 

'±,i 


(m) 


-1 

(4vr ) 2 24F 2 ' 


±,2 


(5.17) 


The replacements given in eq. (5.15) also have to cancel the divergencies in the distribution 
amplitudes for the 5 baryon and the II f and Y f distribution amplitudes, which is the case 
and can be seen as a nontrivial check of our calculation. The higher order divergencies, 
which are contained in our result as a consequence of using IR-regularization [28], have to 
be set to zero by hand. This introduces an unphysical scale dependence in higher order 
terms, which is usually solved by fixing the scale at a typical hadronic value like 1 GeV. 
A variation of this scale within reasonable bounds, say between 0.8 GeV and 1.2 GeV, can 
be used to estimate higher order effects. 

If we neglect higher order contributions, we can rewrite eqs. (5.12) in such a way that 
the complete nonanalytic behaviour is encoded in an overall prefactor: 



(5.18a) 

(5.18b) 
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where 



(5.18c) 

(5.18d) 


- 1 + -AEg • (5- 

From eq. (5.18) it follows directly that at leading one-loop order the complete nonanalytic 
structure is contained in the normalization of the distribution amplitudes, while their shape 
only exhibits the simple dependence on 5m shown in eq. (5.25). Therefore leading finite 
volume effects do only affect the normalization. We want to emphasize that this is only 
true by virtue of our specific choice of DAs. A similar behaviour was found for the meson 
sector (see refs. [29, 30]). The zeroth moments of the given DAs are not independent, due to 
eq. (3.29). In particular all DAs which correspond to operators of certain symmetry classes 
are normalized by the same wave function normalization constants independent of the twist 
of the corresponding amplitude. The zeroth moments define the following normalization 
constants: 



and 


f B = f [dx\$ B i (x 1 ,x 2 ,x 3 ) = ^|§-(l + Asf + )(/* + Sm A/ B ) , 
A? = f [dx]® B 4/5 (x i,x 2 ,x 3 ) = + A gi_)[\l+5m AAf) 

x 2 = f [dx]E b 4 / 5 (xi,x 2 ,x 3 ) = + Ac/f)(A(; + 5m AAf) , 

It = f [dx]Uf(xi,x 2 ,x 3 ) = + A<?n)(/* + dm A/f) , 


It = J [dx]Uf(xi,x 2 ,x 3 ) = yj -^7(1 + Agn)(f* + dm A/f) , 

A t = f [dx] n ^ 5 (xi,x2,x 3 ) = J ^7(1 + A^)(AJ + 5m AAf) . 


(5.20a) 

(5.20b) 

(5.20c) 

(5.20d) 

(5.20e) 

(5.20f) 


For the remaining zeroth moments one finds 


f N = J[dx]Uf (x!,x 2 ,x 3 ) , (5.20g) 

0 = J [dx]$ B 3/6 (x 1 ,x 2 ,x 3 ) = J [dx]E B 4/5 (xi,x 2 ,x 3 ) = J [dx]U 3/6 (xi,x 2 ,x 3 ) , 

(5.20h) 


J [dx]T 4/5 (x 1 ,x 2 ,x 3 ) 



for B * A 
for B - A 


(5.20i) 
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(5.21) 


Due to eq. (5.14), 

A/t = ~ A/ A - l -Af , A/| = ^ A/ a + \ff - A/ w , 

aa| = --aa|--aa| . 

2 2 

In the equations above we have introduced convenient new definitions of / A , Af, Af, /|, 
fj, and A a such that, in the limit of exact SU(3)f symmetry, 

/* = = / S = r = / A = fr = /I , (5.22a) 

A* = Af = Af = Af = Af = A|, (5.22b) 

A 2 = A 2 = A 2 = A 2 = A 2 ■ (5.22c) 

If the reader favors a different definition he or she can easily read off the conversion factor 

from eq. (5.7), noting that additional signs can arise from eq. (5.10) if one uses different 
baryons for the definition of the distribution amplitudes, and that one has to take into 
account additional factors originating from differing definitions of S t , Pi, V), Ai and T* 
(we use the definitions of ref. [ 6 ]). We have performed this matching procedure for the 
constants defined in refs. [5, 6 ] (see appendix C). Note that the constants /|, // and Af 
given above are (at leading one-loop accuracy) completely fixed by f B and Af. However, 
without the knowledge of the SU(3)/ breaking effects one would have to define them as 
additional free normalization constants. /*, A f B , A* and A A f are given by 


f* = f [dx]®l ti (xi,X2,x 3 ) , 

A f B = J [dx] A$ B i (xi,x 2 ,x 3 ) , 

(5.23a) 

A \ = J [dx]$* A/5 (xi,x 2 ,x 3 ) , 

AAf = J [dx]A$> B A/5 (x l ,x 2 ,x 3 ) , 

(5.23b) 

A 2 = J [dx]E* +A/5 (xi,x 2 ,x 3 ) , 

AAf = J [dx] AE b 4 / 5 (xi,x 2 ,x 3 ) , 

(5.23c) 


where, as a consequence of eq. (5.13) (first line) and eq. (3.29) (second line) one has 
Af = -A/ s - Af , AA/ = -AAf - AAf , AAf = -AA 2 S - AA| , 

AAf - -AAf . (5.24) 

The zeroth moments of f 3 ^ and nf/ g (E B 4 ^ 5 and T A ^ 5 ) vanish by construction, since 
they are antisymmetric under exchange of x\ and x 3 (X 2 and x 3 ). One possible approach 
would be to normalize these amplitudes by their first moments. However, our main goal is 
to divide the DAs by normalization constants in such a way that the nonanalytic prefactor 
is canceled. This can be achieved without the definition of additional constants, since all 
prefactors present in eqs. (5.18) also occur in eqs. (5.20). Explicitly, one can consider the 
ratios 



K,i + dm Af - 

->» 

4>* ■ + 5m A4> s ■ 

(5.25a) 

f B 

f* + 5m A f B ’ 

Af 

A* + 5m. AAf ’ 

nf 

+ 5m A<f>ff 

nf 

+ 5m Anf 

(5.25b) 

(5.25c) 

f 

n s/s 

1 

f S/H 

J T 

/* + 5m A f N ’ 

4>; ?; + 5m Anf /S 
~ f*+5m A/| /H 

Af 

A* + 5m AAf ’ 


24 










~B 

1—1 ±,i 

Af 


+ 5m AE ±,i 

A(S + 5m. A\ B 


B 


T 

Af 


J ±B,* 


+ 5m. AT? 


A| + 5m AX B 


(5.25d) 


The idea behind the latter choice is to normalize all DAs with similar behaviour under 
chiral extrapolation (including the ones with vanishing zeroth moment) with the same 
normalization constant containing the complete nonanalytic behaviour. In this way one 
obtains a one-to-one correspondence between a normalization constant and a certain chiral 
behaviour. Note that, following this argument, some of the moments of the leading twist DA 
+ should be normalized with Xf instead of f B . Otherwise the corresponding 
shape parameters do contain chiral logarithms. 


5.3 Example of application 

In this section we will work out explicit expressions for the DAs defined in eqs. (5.7) 
and (5.9) in terms of the shape parameters given in refs. [9-11], where contributions of 
Wandzura-Wilczek type [12] are taken into account explicitly. For brevity we apply the 
approximation advocated in ref. [9], where contributions that can mix with four-particle 
operators are systematically neglected. We use the definitions of said references and we 
define additionally 


'Pnki.Z'l i X 2 j 3 * 3 ) _ Pnk'Pnk(% 3 i X 2 t %l) ■ 


(5.26) 


where p n k = ±1, depending on n and k. This definition is possible since the polynomials 
V n k have definite parity under exchange of x\ and X 3 [9]. We will call the polynomials with 
Pnk = +1 ( Pnk = -1) even (odd). For the DAs we find 
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= Afiff 3 , 
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= Af (iff ^ 
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-B _ x B(^B,WW 4 , ~B,t=5\ 
-±,5~ A 2{- ± ,5 +- ± ,5 )’ 


(5.27a) 

(5.27b) 

(5.27c) 

(5.27d) 

(5.27e) 

(5.27f) 


where all chiral logarithms are contained in the prefactors. Analogous expressions for the 
II and T DAs will be given below in eq. (5.34). Genuine twist 5 contributions (<3?^g -5 , 
will be neglected in this approximation. Also twist 6 DAs are neglected; one could 
in principle take into account Wandzura-Wilczek contributions to the twist 6 DAs, but 
the corresponding expressions are not known yet. The shape of the DAs is given by the 
genuine twist 3 and twist 4 contributions 


3 (xi, X 2 ,X 3 ) = 120 xiX 2 X 3 Y l l PnkPnk(xi,X 2 ,X 3 ) , (5.28a) 

n,k<n 
Pnk 1 

® B £ = 3 ( X 1, X 2, X 3) = 120xiX2X 3 Y l l PnkP n k(xi,X2,X3) , (5.28b) 

n,k<n 
Pnk = ~1 


- 25 - 





$ 


4- 


B,t= 4/ 

+,4 l 

B,t=4 

—,4 

B,t=4 




(xi,x 2 ,x 3 ) = 24xix 2 (?7,fo + 2(2 - 5x 2 )??f 0 
x 2 ,x 3 ) = 24x 2 x 3 ^f 0 -jjo 


— B ’4 4 (x l5 x 2 ,x 3 ) = 24x 2 x 3 ^-^(x2-x 3 )^fo + ••• j , 


V 11 + • • • ) , 

(5.28c) 

+ ...), 

(5.28d) 


(5.28e) 

4 

(5.28f) 


and the Wandzura-Wilczek contributions (see refs. [9-11]) 

240^ B fc ( _ 9 
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X XiX 2 X 3 V n k(xi,X 2 ,X 3 ) , 

^+’^ W 4 (xi,x 2 ,x 3 ) = 4x 3 (5(x 4 + 2x 2 x 3 -Xg)?/^ + ...) , 
^f W4 (®i,x 2 ,x 3 ) = 4 x 3 (1 - x 2 )(2rjQQ + 3(1 - 5x 2 )??f 0 + ...) , 


= 4x 1 (1 + xi)^ 0 - ?(4-4xi +x?-5x?Kf 0 + 


-+,5 


5^(11,12,13) = -12 xi(x 2 - x 3 )Cto + W “( 5 - x \ + 5 xf)(x 2 - x 3 )^f 0 + ... 


(5.29a) 

(5.29b) 

(5.29c) 

(5.29d) 

(5.29e) 

(5.29f) 

(5.29g) 

(5.29h) 


where the summation over n starts from 0 and, generally, goes to infinity, but is truncated 
at n = 2 in the approximation of ref. [9]. 1 Note that our separation into “+” and ” 
amplitudes at leading twist level corresponds to a separation of even and odd polynomials. 
The normalization constants are still defined such that i)q 0 = <p B 0 = £ B 0 = 1, which are only 
kept for a cleaner notation. Note also that the introduction of (p B , and fj B k only amounts to 
a redefinition of the shape parameters occurring in anc [ W3 by a factor 

of f B /Xf and the ones occurring in 4>^ 4 -4 and <&b,WWa ^y a f ac tor of Xf /f B with respect 

1 We do not take into account possible quark mass corrections to eq. (5.29) and eq. (5.38) below (compare 
e.g. refs. [31, 32] where such computations have been performed for vector-meson and pseudoscalar-meson 
DAs), since they can (by definition) be absorbed into the genuine higher twist terms in eq. (5.27) and (5.34). 
Let us note in passing that our general result does not rely on the separation of Wandzura-Wilczek and 
genuine higher twist terms at all, since the calculation within chiral perturbation theory does not distinguish 
between these contributions. 
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to ref. [9] (the corresponding anomalous dimensions have to be adjusted accordingly): 


f B \ B 

~B J B ~B ^1 B /r q r\\ 

Vnk ~ ( P‘nk ? Vnk ~ j q Vnk * (5.30) 

The dependence of the shape parameters on the quark mass splitting takes the following 
form 
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if Pnk = -1 , (5.31a) 

B VlO + Ar ll0 

7/10 A* + 5m AAf ’ 

(5.31b) 


(5.31c) 


which corresponds directly to eq. (5.25), while the dependence of the normalization con¬ 
stants is given in eq. (5.20). The parameters describing SU(3)j symmetry breaking are 
restricted by eq. (5.13) such that 


Ax nk = ~ Ax nk ~ Ax nk > for ^ 6 W, <P,V, ■ (5-32) 


For the original twist 3 and 4 DAs given in ref. [6] (see also eq. (5.6)) the new choice of 
normalization yields 
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where, as discussed above, the normalization of the odd moments of the leading twist 
amplitude with Af (instead of f N ) appropriately reflects their chiral behaviour. Note that 
this is consistent with an earlier two-flavor BChPT calculation, where it was found that 
the odd first and second moments of the leading twist amplitude have the same chiral 
logarithms as Af (see appendix of ref. [13]). 
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For a description of the complete baryon octet one also needs the II and T DAs defined 
in eq. (5.9), which are relevant for the hyperons. These are completely fixed by the 4> ± and 
H DAs. Consequently, the following equations do not contain any additional parameters: 
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(5.34a) 
(5.34b) 
(5.34c) 
(5.34d) 
(5.34e) 
+ Tf’* =5 ) , (5.34f) 
(5.34g) 


where the genuine twist 5 contributions 11^’* 5 and Y,? 4 5 will be neglected as above. The 
genuine twist 3 and twist 4 contributions are 
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T 


) 

)■ 


The shape parameters are fixed: 


7T, 


s/s Att^ 

nfe /* + A/^ /h 

~s/s = f?ii + <5m Aff/ a 
11 /* + <5m A/^ /H 


vr: 


A 


nA: 


nA: 


A 4 + dm AA^, 


A = 7?io + 5m ACiq 
^ 10 A* + <5m AA^ ’ 


v?o = 


£i 0 + <5m An 


10 


A* + 5m. AAf ’ 


(5.35a) 

(5.35b) 

(5.35c) 

(5.35d) 

(5.35e) 

(5.35f) 

(5.36a) 

(5.36b) 

(5.36c) 


Y' I'-' A 

where A f T ~ and AA^ are defined in eq. (5.21). The parameters describing SU(3)/ sym¬ 
metry breaking can be determined by eqs. (5.13) and (5.14): 


A?r nfc - ~2 ~~ 2^-Vnk. 


A^ = -±A^~A¥&, 


- 28 - 







A ^nk = \ A Vnk + \ A ^lk ~ A( Pnk , 

A Cn = ~ l A Vu , 

A Cn = 2 A Vii + 2 A Vn ~ A Vn > 

At, io = ~2 A €io ~ 9 A £lO > 

Au ib = ^ A £io + 9 A £io ~ A £io ■ 

The Wandzura-Wilczek contributions take the form 

E/S 


A Cu) = -l A Vio - l Ar ll 0 


Au 10 “ 9 A ?10 9 A £l 0 ) 


n 


s/s ,ww 3 


240vr 

(xi,X2,x 3 ) = -£———— 
n %in + 2 ){n + 2 ,) 


n 


A ,WW 3 


(X 1 ,X 2 ,X 3 ) = ~Y, 


24074 


Pnk~ 1 


(n + 2 - )xiX2.'C3'P nfc (xi,X2,.T 3 ) 

(' " 8 


, o\/ 0 s,n + 2- -—]xix 2 x 3 V n k(xi,X 2 ,x 3 ) 

k<n\ n + 2)( n + 3) \ 3 x 3 ; 


n 


E/S ,M/W 3 


(xi,x 2 ,x 3 ) = )T 


240vr 


S/2 

nk 


n,k<n 
Pnk = ~*~ 1 


(re + 2)(re + 3) 

x x 1 x 2 x 3 V nk (xi,x 2 ,x 3 ) , 


7 3 \ 

l 3 \ . ,7 

K" + 2 '8F) 



H- 3 ’ WW 3 ( X l,X 2 , X3 ) = X! 

n,k<n 
Pnk = ~1 


240vt^ 


(re + 2 )(re + 3) 

X XlX 2 X 3 P nk (x 1 ,X 2 ,X 3 ) , 


i n+2 -^)(" +1 -^)- 


(re + 2) 2 


n s/=,u w *( XuX2}X3 ) = 4x 3 (5(xi + 2x 2 x 3 - x 3 )Cf/“ + ...) , 
H 3 ' WWa ( X i-, X 2 -, X z) = 4x 3 (1 - x 2 )(2Cto + 3(1 - 5x 2 Km + • • • ) - 
T?/ _J1 W4 (xi,x 2 ,x 3 ) = 4xi(1 + xi)uq 0 / “ - y(4-4x! + x\ - 5xi)re E/ ““ 


tA,wW4(xi ) x 2) x 3 ) = -12xi(x 2 -x 3 )uqo + y(5 ~x± + 5xf )(x 2 - x 3 )uf 0 + ... . 


10 + • • • 5 

A 


(5.37a) 

(5.37b) 

(5.37c) 

(5.38a) 

(5.38b) 

(5.38c) 

. (5.38d) 

(5.38e) 

(5.38f) 

(5.38g) 

(5.38h) 


To conclude this section we want to point out the merits of our calculation. First of all, 
we found that the behaviour under chiral extrapolation of a certain moment correlates to 
its parity in the sense of eq. (5.26). Therefore it is advantageous to normalize the odd 
moments of the leading twist DA with \f instead of f B . Quantitatively more important, 
however, is the significant reduction of parameters: we find that (within the approximation 
used above) we only need 43 parameters to describe the complete set of baryon octet three- 
quark DAs (including their dependence on the quark mass splitting). In contrast, an ad hoc 
linear extrapolation without the knowledge of SU(3)y symmetry breaking would require 
72 parameters for the given setup, since one can not make use of eqs. (5.21), (5.24), (5.32) 
and (5.37). 


5.4 Dependence on the mean quark mass 

The distribution amplitudes <f>* % and H* t have a nontrivial dependence on the mean quark 
mass fhq. This is not really interesting from a phenomenological point of view, since the 
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number of independent distribution amplitudes can not be further reduced compared to 
eq. (5.18), even if one expands everything around the chiral limit. However, the dependence 
is of importance for the analysis of lattice data if one wants to include data points from 
simulations with unphysical mean quark mass. The mass dependence reads 


^ = ^(l + ^'* +5 J ± ) + mA<hV , 
SV = S; ): (l + iS'* + fe) + mAHV, 

where 

l 2 Bofh q 


(5.39a) 

(5.39b) 


(5.40) 


g^ ± , gl and T,'* are functions of the mean quark mass that can be taken from appendix A. 
The divergencies occurring at linear order in the mean quark mass can be canceled via the 
following introduction of counterterms 


A<b 


mfc * 


±,Z 


24 F? 


f^AT + A^r^), ASA 


rn b C E t \ —+ ,ren. 

24F* ^ ±4 


(A*) > ( 5 - 41 ) 


where L contains the divergence (see appendix A) and the coefficients are 


4± = -(6(5T> 2 + 9 F 2 ) + 13 ± 6) , 
3 


c| = -(6(5T> 2 + 9F 2 ) + 9) . (5.42) 

3 


This leads to the following scale dependence in the renormalized amplitudes: 


9 A i * ren. / \ 1 ^<I>± ,*.0 9 A ^* ren. / . 

V^ ; ‘ ^ ‘ to- 


-1 mf cl 


b 


(4vr) 2 24F 2 




(5.43) 


The divergencies occurring together with higher orders of the quark masses have to be 
canceled by hand as discussed in section 5.2. If one takes eq. (5.39) and plugs it into 
eq. (5.18) one finds (up to terms of higher order) 


<i = \/^(l+4 ± + A 5 £ ± )(<]?A + m A$l >i + 5m A<J?y , (5.44a) 

-±,i = \/^( 1+ 9s + A ai)(K,i + ™ ASA + <5m ASf^) , (5.44b) 

H? = y/z^( 1 + 9L b + + fn A+ dm AHf) , (5.44c) 

Tf = ^(l+ 5 3 + A 5 |)(s: Bii + m A~l Bti + 5m AT,f) , (5.44d) 

where 

Jz^=\ + hl'* + ^' B . (5.45) 


Starting from this point everything can be worked out analogously to the case of fixed 
mean quark mass. 
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6 Summary 

In this work we have presented the first analysis of baryon octet light-cone DAs in the 
framework of three-flavor BChPT. At next-to-leading order accuracy in the chiral counting 
scheme, we obtain the leading quark mass dependence and (automatically) the leading 
SU(3)/ breaking effects. Describing the baryon octet simultaneously we are able to unify 
and systemize the efforts made in refs. [5-8, 33]. 

An important insight to be gained from our results is of qualitative nature: in the 
chiral odd sector the chiral behaviour (i.e. the contained chiral logarithms) of a specific 
moment does not depend on its twist, but on whether it contributes to the or $ B i 
amplitudes (see eq. (5.7)). Those contributing to the “+” (“-”) amplitudes have the same 
chiral logarithms as f B (Xf ). Therefore the odd moments of the leading twist DA behave 
like Af instead of (as one might have expected) f B . This result is consistent with an earlier 
two-flavor calculation, where it was found that the odd first and second moments of the 
leading twist DA have the same chiral logarithms as X B ' (see appendix of ref. [13]). 

In section 5 we provide a set of DAs that parametrize the complete baryon octet 
(including the A baryon) in a minimal way and do not mix under chiral extrapolation. 
Eqs. (5.18) and (5.44) are our main results. They describe the quark mass dependence of 
the baryon octet DAs (including all higher twist amplitudes) in a very compact manner. 
Eq. (5.20) contains explicit extrapolation formulas for the wave function normalization 
constants, while the dependence on the quark mass splitting of the shape parameters, 
which describe contributions of higher conformal spin, is shown in Eqs. (5.31) and (5.36). 
The results will be of particular importance for the interpretation and extrapolation of 
forthcoming lattice QCD data, due to the significant decrease in number of parameters 
(compare section 5.3). For the same reason our results are relevant for QCD sum rule 
analyses and for model building. 
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A Loop contributions 

The functions g^ A and A g B A are given by 

6 F?gl + = -19 H^) + 2(5 D + 6F)R 2 (m* m ) , 

§F*g%_ = -7Hi - 10DH 2 (m^) , 

6 FlgX = - 18FH 2 (nO , 

24F* Ag$ + = -57 AHi(?n jr ) - 18 AHi(m^) - AHi(m r? ) + 30(D + F) AH 2 (m n ) 

+ 12(D + F ) AH 2 {rriK) + (-2 D + 6 F) AH 2 (m^) , 

24F* A</$ + = -12 AHi(m 7r ) - 60 AHi(mA') - 4AHi(m ?? ) + 24D AH 2 (m 7r ) 

+ 24 (D + 2 F) AH 2 (mft-) - 8 D AH 2 (m ri ) , 
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24F 2 Ag$ + = -9 AHi(m„-) - 66 AHi(mx) - AHi(m r? ) + 18 {-D + F ) AH 2 (m n ) 

+ (6014 + 36T 1 ) AH 2 (m^) - 2(14 + 3 F) AH 2 (m v ) , 

2AF+ Ag$ + = -36 AHi(?n 7r ) - 36 AHi(mA') - 4AHi (m v ) + 2414 AH 2 (?n 7r ) 

+ 8 {D + 6 F) AH 2 (mx) + 8 D AH 2 (m r) ) , 

24F 2 Ag$_ = -9 AHi(m 7r ) - 18 AHi(toa') - AHi (m v ) - 18(14 + F) AH 2 (m n ) 

+ (-20 D + 12F) AH2(?7 ir) + {-2D + 6 F) AH 2 (m^) , 

24F 2 Ag$_ = -12 AHi(m 7r ) - 12 AHi(?ba') - 4AHi(m^) - 814 AH^m^) 

- 2414 AH 2 (mi<-) - 814 AH 2 (m ?7 ) , 

24T l2 Agf_ = -9 AHi(m 7r ) - 18 AHi(toa) - AHi(?n ?? ) + 18 (-D + F) AH^m^) 

- 4(5 D + 3 F) AH 2 (mA') - 2 (D + 3 F) AH 2 {m v ) , 

24F 2 Ag$_ = -36 AHi(m 7r ) + 12AHi(?nA-) - 4AHi {m v ) - 72D AH 2 (m 7r ) 

+ 2414 AH- 2 {rriK) + 814 AH 2 ( 771 , 7 ) , 

A g% = A g% + , 

24F+Ag B = -24 AHi (m n ) - 36 AHi(mA') - 16 AHi(m^) + 48-F AH 2 (m„-) 

+ 24 D AH 2 {mx) + 16 D AH 2 (m^) , 

24F 2 Agg = -9 AHi (m n ) - 42 AHi(toa) - 25 AHi(m^) + 18(14 - F) AH 2 (m,r) 

+ 12 (D + 3 F) AH 2 (mA') + 10 (D + 3 F) AH. 2 {m ri ) , 

24F 2 A^n = _ 12 AHi {thr) - 16 AHi(m 7? ) - 24 D AH 2 (mA') - 16 D AH 2 ( 771 , 7 ) , 

24F 2 AgS = -9 AHi(m 7r ) - 18 AHi(mA') - 9 AHi(m r) ) - 18(14 + F) AH 2 {m^) 

+ 12(14 - 3 F) AH 2 (mA') + 6(14 - 3 F) AH 2 (m, 7 ) , 

24T’ 2 Ag= = -24AHi(?n7r) - 12 AHi(mA') - 48F AH 2 (m vr ) - 24 F AH 2 (m^) , 
24F+Ag§^ = -9 AHi(m 7r ) - 18 AHi(mA') - 9 AHi(m r) ) + 18(14 - F) AH 2 {m^) 

- 12(14 + 3 F) AH 2 (toa) - 6 {D + 3 F) AH 2 ( 771 , 7 ) , 

2AF?Ag£ = -36 AHi(mjr) - 72F AH 2 (m A -) . (A.l) 


The Z-factor contributions are given by 


S'* = ^(5H 2 + 9F 2 )H 3 (m^) , 


AT,' b = 3 g B ,n AH 3 (m 7r ) + 4 g B j< AH 3 {m K ) + 9b, AH 3 ( 771 ,,) , (A.2) 


where the coefficients gB,M are defined in eq. (4.7). The auxiliary functions AH k are defined 
as 


AH k (m) = H k (m) - H k (m^) 


(A.3) 


with 



(A.4a) 
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H 3 (m) 
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32-K 2 mt 2 8ir 2 m) 
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rrr 


4 m f 


■arccos 


(m \ 

"M/ ’ 


/ 2m 2 \ [ 1 / m 2 \ 

l 1_ 3^jr + 32^ l0g ('^'j. 


m 


32F 2 vr 2 


3m d 


I 1 '#) 


32F 2 m?ir 2 




: arccos 


4mf 


\ 2 K / 


(A.4b) 


(A.4c) 


L contains the divergence and the finite constants typical for the modified minimal sub¬ 
traction scheme in 4 - e dimensions: 


L E 


(e + ^( 1 + log ( 47r ) - ^)) ' 


(47t) 2 


(A.5) 


Note that we have shown that the divergencies of leading one-loop order can be canceled. 
For practical purposes one can therefore set L to zero everywhere if one simultaneously 
replaces the corresponding DAs by the renormalized ones (compare section 5). Within our 
level of accuracy it is legitimate to replace m£ and F 1 * by their values at the symmetric 
point, where m q = fhq hys . 


B Handbook of distribution amplitudes 

In this section we express the the 24 standard DAs occurring in the general decomposition 
derived in ref. [6] (S B , P B , V B , Af, T B ) in terms of the DAs defined in section 5. The 
equations given below follow directly from the definition of the DAs in eqs. (5.7) and (5.9) 
together with the symmetry properties of the standard DAs under exchange of the first 
and the second variable given in eq. (4.21). For the twist 3 and twist 6 amplitudes one 
finds 


V, 
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r/6- 2 
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\ C B C ~B ) 
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l C B C B ) 


H B I§{ X 1,X3,X2) 


(-1) 


B " 


(B.l) 


-'B 


where the DAs on the l.h.s. are functions of (aq, a; 2 , aq). The twist 4 and twist 5 amplitudes 
read 
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(B.2) 


C Matching to other definitions in the literature 

Since we use the same definitions as ref. [6] it is no surprise that 

f N = f N ([ 6]), Af = A 1 ([6]), A? = A 2 ([6]). 


(C.l) 


We can also match some of our constants to the leading twist normalization constants given 
in ref. [5]. Note that for the E and S a relative minus sign originates from the fact that 
ref. [5] uses E + and H _ for the definition, while our choice is XT and E° in order to have 
the same sign as for the proton. 


f = /jv([5]) , 
f S = -/e([B]) 
f = -/e([5]) 
^2 


It ~ ~/s([5]) , 

/! = -/!([ 5]) , 


(C.2a) 

(C.2b) 

(C.2c) 


/ =\/3/a([5]) 


J [dx]x 1 $ A 3 (x 1 ,x 2 ,x 3 ) = n/6/a ( [5] ) • (C.2d) 
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Due to some misprints, obvious errors and inconsistencies within refs. [7, 8 , 33] we are not 
able to give reliable matching formulas for their definitions. 


D Some construction details 


In the first part of this section we will describe why we can trade covariant derivatives acting 
on the baryon field for normal derivatives acting on the complete current. This choice is 
very convenient since the external derivatives (in contrast to the covariant derivatives acting 
on the baryon field) do not lead to additional loop momenta in the integrals. To show that 
this formulation only differs in higher order contributions we use the identities 


e 

£ 


abc a'b'c' 

— ^aa'^bb'^cc'^ 


abc = ^) aa ,{u^ w {v)) cd e a ' b ' c ' , 

0 = ( (<9/^)aa' ^bb'^cc' "*■ Uaa' ( ) bh' U-cc' ^ aa' ^bb' (^jl ^ ) cc J ) - ■ 

0 = )aa'(^ )bf/)cc' )feb' (u )cc' )aa'(^ 


(D.la) 

(D.lb) 

(D.lc) 


(D.ld) 


which follow from det(u) = 1. From these one obtains 

(- D, L B) aa ,e a ' bc = ({d^B) aa ,5 w 5 cc , + {T^B) aa ,5 w 5 cc , - (BT^) aa ,5 w 5 cc ,)e a ' b ' c ' 

= {(d,B) aa' fibb'^cc' “I" (r fiB^aa'^bb'^cc' B) aa ' (T^ 5 CC ' (D.2) 

+ {B) aa ,5 w {T il ) cc ,)e a ' b ' c ' . 

Additionally we need 

d^ux = u x (uxd^u x ) = rix(r M - (-1 )x^u^j = u x T^ , (D.3a) 

9/jXm = D^Xm - [T M ,Xjvf] = -[T M ,Xjvf] , for Xm e {u v ,x±} ■ (D.3b) 


Putting everything together we find for a general structure with arbitrary mesonic building 
blocks Xi,X 2l X 3 € {u u , x±} 

( u x XiD fi B) aa '(uYX 2 )bb'(uzXs) cc '£ a b c 

= ((uxX^B) aa' (uYX 2 )bb'(uzXs)cd ('U'XX\l? aa ' (llyX2(j-^zX^)cc' 

+ (u x X\B) aa / ( uyX-2 r /( ) w ( uzX :i ) cc i + ( u X XiB) aa '(uYX2)bb'(uzXsT fl ) cc ')£ a b c 
= d,j((ii x X] B) aa i(uyX 2 )bb'( u zX;$) cc >)e a h c . (D.4) 

In the last step we have used 


uxXiTj, = u x d^X 1 + u x T^X 1 = u x dfj,Xi + {d^u x )Xi = d^{u x Xi) , (D.5) 

and the same for uyX 2 and uzX 3 . 

In the following we will argue that structures involving baryon and meson fields at 
different positions can be dropped. We can choose the structure containing the baryon to 


- 35 - 


be situated at x, while we call the second position y such that we can write schematically 
B(x,y ) = f(x)g(y), where g only contains mesonic building blocks. Every derivative acting 
on g therefore has to be counted as first order in the chiral power counting. It follows 
trivially that 

f(x)g(y) = f(x)(g(x) + {x-y)- dg{x) + ...) = f(x)g(x) . (D.6) 


E Matching relations 

In this section we provide the result of the matching described in section 4.3, which is needed 
in intermediate steps of our calculation (in practical applications one can always use the 
readily evaluated expressions shown in appendix B). For the twist projected amplitudes 
introduced in ref. [6] one finds 
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bca . , .7 l, cab , r,i2,cab , 7 4, cab , 

,od^ X2 ' X3,X1 ’ + B,odd'- X3,xl,X2 ' ~ ltl B,odd( X3 ’ xl ’ X2 > + 4 /l B,odd( X3,xl ’ X2 ' 


2 , bca 


4, bca 


7, bca 


2 K 


B 


+ 2/l^“ dd (LC 3 ,xi,a- 2 ) - 2h£“ dd (x 3 ,xi,x 2 ) , 

4 ,bca Qi7,bca , x A ^8,bca 


8,cab 


V d = 8h5^ dd (x 2 ,x 3 ,xi) + 8 h^ dd (x 2 ,x 3 ,xi) - 4 h^“ d (x 2 ,x 3 ,xi) - 16h^ dd (x 2 ,x 3 ,xi) 

+ 8/l^“ dd (x 3 ,Xl,X 2 ) + 8h^ dd (x 3 ,xi,x 2 ) - 4/l^“ dd (x 3 ,Xl,X 2 ) - 16/l^“ dd (x 3 ,Xl,X 2 ) , 

aB a 1 8 , cab . aiS, bca . 

A l = 4/l Bi0dd (x 3 ,*i,* 2 ) - 4/l B,odd( X2 ’ X3,Xl ' ’ 

A 2 = 2h B ^ dd {x2,X 3 ,X 1 ) + 2h 7 B ^ dd (x2,X3,Xl) - 2h 8 B b ^ dd (x 2 ,X 3 ,Xl) 


+ 2 h 


2,cab 
B 


cau , x 0 7 7 ,cab . , 0 7 o ,cau , » 

oddl® 3 ^ 1 ’ 212 ! " Zh B,odd^ X3 ’ xl ' X2 > + Ztl B,odd( X3 ’ xl ’ X2 > 

A B = 2h 2 B b Q dA (x2,X 3 ,Xl) - 2h 7 B bC odd {x2,X 3 ,X 1 ) + 2/lg 6 (( dd (x2,X 3 ,Xl) 


8, cab 


+ 2 h 


2,cab 
B 


cao . . . 0 7 7 ,cab . . 

,odd' X3,Xl,X2 ' + 2n B,odd( X3 ’ Xl ’ X2 l 


c\i&,cab . . 

Ztl B,odd ( ' X3 ’ Xl ’ X2 > ’ 
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a b 

Sl 4 


1 B 


A* = 


t b 

rjiB 
1 2 


~^ 1 B,odd ( ' X2 ’ X3 ’ x A + ^B,odd( X2 ’ X3 ’ Xl>) ~ ^Bfidd^ 2 ' X3 ’ X A ~ ^Bfid d( x 2,x 3 ,xi) 

+ ^BMd^ X2 ' X3 ' Xl ) ~ Atlj^^ A {X3,Xl,X2) + 2h 2 g C ^ d (x 3 ,x 1 ,x 2 ) + 4/l^“j d (X3,Xl,X 2 ) 

+ ^B^oddfx3 ,xi,x 2 ) - 2h 8 g C ^ d (x 3 ,xi,x 2 ) , 

~^ 1 B b ,odd ( ' X2 ’ X3 ’ x A + ^ 2 J > ^ dd (x 2 ,X3,xi) + 4:h^ dd (x2,X3,Xl) + 2h 7 ^ dd {x 2 ,X3,xi) 

- ^h 8 J > ^ )dd (x2, x 3, x i) - 4/l^“j d (x 3 ,a:i,a:2) + 2/l£“ dd (x 3 ,xi,X2) - 4/j.^“j d (® 3 ,xi,X2) 

- 2^^Q dd (x 3 >Xl|X2) + 2/l^“ dd (x 3 ,Xl,X 2 ) , 

^h^ dd (x 2 ,x 3 ,x 1 ) + 8/l^ dd (x 2 ,X 3 ,Xl) - 4/j.^ dd (X2,X 3 ,®l) - 16/l^ dd (x 2 ,X 3 ,Xl) 

- 8/l^“ dd (x 3 ,a:i,X2) - 8/l^“ dd (x 3 ,®i,x 2 ) + 4/l^“ dd (x3,®i,® 2 ) + 16/l^“ dd (x 3 ,xi,x 2 ) , 

: 4 ^'B^odd ( X1 ' X2 ,X3 ) > 

: — 1 6/l^^^ en (xi,X2,X 3 ) , 


T. 


■B 


T, 


B 


rpB _ 

1 5 ~ 

nB 


~^B^ven^ Xl ’ X2 ’ x A + ^B^eve n ( xi >®2,x 3 ) ~ 2/l^“ fe dd (xi,X2,x 3 ) + 2/!^“ 6 dd (o;i,a:2,x 3 ) , 

_4/l B a even(* 1 ’ a:2 ’ a;3 ) “ ^fCeven^ 1 -* 2 -* 3 ) “ ^“even ( x UX2,x 3 ) ~ 2h 6 ^ e J X1 ,X2,x 3 ) 
~ 4^ a o b dd (xi,x 2 ,x 3 ) - 2/l^“ b dd (xi,x 2 ,x 3 ) + 2/l^“ fc dd (a;i,®2,x 3 ) , 


q 7 4,d6c . . -1 r* l&jdbc . . 

^B,even^ Xl,X2 ' X A ~ ibh B,even^ X2 ' x A , 


f 5,afrc 


= - 8 /i 


4,abc 
B, odd 


T- 


B 


Tc 


B 


(xi,X2,x 3 ) - 8/?^ a o fe dd (xi,a;2,X3) + 4 ^B^dd( xl > X2 > X3 ) + 16/i£“o dd (xi,X2,x 3 ) , 

0 i5,abc . . 0 7 6,a6c . . 0 7 7,abc . . . 0 7 8,a6c . . 

OB 3 even \ x l^ x 2, x 3) - Zn B,eveiA Xl ’ X2 ’ X3 > ~ Zr> ‘B, odd^ 1 ’* 2 ’ 3 ' 3 ) + ^B.odd^ 1 ’ 352 ’® 3 ) ’ 

^bJLi^ 1 ’* 2 ’*®) + ^^n^ 1 ’* 2 ’* 3 ) + ^B^L^ 1 ’* 2 ’* 3 ) + ^“even^ 1 ’* 2 ’* 3 ) 


■4 h 


4 ,abc , r-r 7,a6c , . r> 18 ,afec , 

B,odd( X1 ’ X2,X A ~ ^B.odd^ 1 ’® 2 ’^ 3 ) + ‘^B.oddl 351 ' 12 ’* 3 ' ) 


(E.l) 


where the DAs on the l.h.s. are functions of (xi,X 2 ,X 3 ). The functions on the r.h.s. are 
given in eq. (4.17). For the flavor indices a,b,c on the r.h.s. one has to insert the flavors 
of the operators for which the l.h.s. is defined. A standard choice is p = uud, n = ddu, 
E + = uus, S° = uds, T,~ = dds, S° = ssu, E~ = ssd, A = uds, where the order of the flavors 
is relevant for the symmetry properties of the DAs. 
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